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We discuss exclusive production of lepton-antilepton, quark-antiquark, tvtt and p°p° pairs in ultraperipheral, 
ultrarelativistic heavy-ion collisions. 

The cross sections for exclusive production of pairs of particles is calculated in Equivalent Photon Approx- 
imation (EPA). Realistic (Fourier transform of charge density) charge form factors of nuclei are used and the 
corresponding results are compared with the cross sections calculated with monopole form factor used in the 
literature. Absorption effects are discussed and quantified. The cross sections obtained with realistic form factors 
are significantly smaller than those obtained with the monopole form factors. 

The cross section for exclusive p + pT production in nucleus - nucleus collisions are calculated and some dif- 
ferential distributions are shown. The effect of absorption is bigger for large muon rapidities and/or large muon 
transverse momenta. We present predictions for LHC. 

We calculate cross section for exclusive production of tt + -k~ and 7r°7r° pairs. The elementary process 77 — > tttv 
is discussed in detail. We concentrate on high-pt processes. We consider pQCD Brodsky-Lepage processes or 
alternatively hand-bag mechanism. The nuclear cross section is calculated within b-space EPA for RHIC and 
LHC. 

Similar analysis is performed for p°p° production, where the elementary cross section is less known. Our analysis 
includes a close-to-threshold enhancement of the cross section. The cross section for the low-energy phenomenon 
is parametrized and the high-energy cross section is calculated in a simple Regge model. Predictions for heavy 
ion collisions are presented. 

The cross section for exclusive heavy quark and heavy antiquark pair (QQ) production in heavy ion collisions is 
calculated for the LHC energy ^/sjvjv = 5.5 TeV. Here we consider only processes with photon-photon interactions 
and omit diffractive contributions. We include both QQ, QQg and QQqq final states as well as photon single- 
resolved components. The different components give contributions of the same order of magnitude to the nuclear 
cross section. The cross sections found here are smaller than those for the diffractive photon-pomeron mechanism 
and larger than diffractive pomeron-pomeron discussed in the literature. 



1. Introduction 

It was shown in several review articles [T] that 
the ultrarelativistic collisions of heavy ions pro- 
vide a nice opportunity to study photon-photon 
collisions. This is due to the enhancement caused 
by the large charge of the colliding ions. Paramet- 
rically the cross section is proportional to Z\Z\ 
which is a huge number. It was discussed recently 
that the inclusion of realistic charge distributions 
in nuclei lowers the cross section compared to the 
naive predictions. Recently we have studied the 



production of p°p° pairs [2J, of fi + pT pairs [3J, of 
heavy-quark heavy-antiquark pairs [1] as well as 
of DD meson pairs [5] . 

Here we shall briefly summarize the recent 
works. 



2. Formalism 

2.1. Equivalent Photon Approximation 

The equivalent photon approximation is a stan- 
dard semi-classical alternative to the Feynman 
rules for calculating cross sections of electromag- 
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Figure 1. A schematic picture of the collision and the 
quantities used in the impact parameter calculation. 



netic interactions This picture is illustrated 
in Fig. Q] where one can see a fast moving nucleus 
with the charge Ze. Due to the coherent action 
of all protons in the nucleus, the electromagnetic 
field surrounding (the dashed lines are lines of 
electric force for a particles in motion) the ions is 
very strong. This field can be viewed as a cloud of 
virtual photons. In the collision of two ions, these 
quasireal photons can collide with each other and 
with the other nucleus. The strong electromag- 
netic field is a source of photons that can induce 
electromagnetic reactions on the second ion. We 
consider very peripheral collisions i.e. we assume 
that the distance between nuclei is bigger than 
the sum of radii of the two nuclei. Fig. Q] explains 
also the quantities used in the impact parameter 
calculation. We can see a view in the plane per- 
pendicular to the direction of motion of the two 
ions. In order to calculate the cross section of a 
process it is convenient to introduce the following 
kinematic variables: 

• x = lu/Ea, where uj energy of the photon 



and the energy of the nucleus 

• Ea = lArriproton = jMa where Ma is the 
mass of the nucleus and Ea is the energy of 
the nucleus 

The total cross section can be calculated by the 
convolution: 

a (AA -> cic 2 AA; s A a) = 



a (77 -> cic 2 ; W ni = ^xix 2 s A a) 



An 



77 (a;i,X2,b) . 



(1) 



The luminosity function can be expressed in 
term of flux factors of photons prescribed to each 
of the nucleus. 



dn 77 (wi,w 2 ,b) 



S 2 abs (b) 



d 2 b 1 7V(wi,b 1 )d 2 b 2 iV(a; 2 ,b 2 ) 



dwi du>2 

U>1 LO2 



(2) 



The presence of the absorption factor S^ bs (b) 
assures that we consider only peripheral colli- 
sions, when the nuclei do not undergo nuclear 
breakup. In the first approximation this can be 
taken into account as: 



6(b-2R A ) = 6(\b 1 -b 2 \-2R A ) . (3) 



S 2 abs (b) 

Thus in the present case, we concentrate on 
processes with final nuclei in the ground state. 
The electric field force can be expressed through 
the charge form factor of the nucleus [3] . 

The total cross section for the AA — > C\C 2 AA 
process can be factorized into an equivalent pho- 
tons spectra and the 77 — > cic 2 subprocess cross 
section as 

a {AA -> c\c<iAA; s A a) = 

J a (77 cic 2 ; V4wiw 2 ) 9 (|b>i - b 2 | - 2R A ) 

iV( Wl ,b 1 )7V( W2 ,b 2 )d 2 b 1 d 2 b 2 ^i^ 2 - . (4) 

We introduce the invariant mass of the 77 sys- 
tem: W77 = Additionally, we define 
Y = I (y Cl +Vc 2 ) rapidity of the outgoing cic 2 
system which is produced in the photon-photon 
collision. Making the following transformations: 



77 P Y 



it: 



w 2 = 



77 -y 



(5) 
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dwi dw 2 



-dW 77 dF 



dwidw 2 -> dW 77 dY 



(6) 



(7) 



formula (HJ) can be written in an equivalent way 

as 



cr (AA -> cic 2 AA; s AA ) = 
<r(77^ cic 2 ; W 77 )6>(|bi 



2R A ) 



JV(wi,bi)JV(wa,ba) 



xd 2 b!d 2 b 2 



-dVT 77 dF 



(8) 



77 



Finally the cross section can be expressed as the 
five-fold integral: 

a ( 77 -> W / 77 ) (|bi - b 2 | - 2fl A ) 

xA^(o;i,bi)iV(a; 2 ,b 2 ) 



27r& m d6 m dfr-r d&, 



77 dJT 77 dF : 



(9) 



where b x = (b lx + b 2x )/2, b y = (b ly + b 2y )/2 and 
b m = b\ — b 2 have been introduced. The formula 
above is used to calculate the total cross section 
for the AA — > AAc\c 2 reaction as well as distri- 



butions in b 



M ClC2 and Y(cic 2 ) 



Different forms of form factors are used in the 
literature. We compare the equivalent photon 
spectra for realistic charge distribution and for 
the case of monopole form factor. 

2.2. Charge form factor of nuclei 

The charge distribution in nuclei is usually ob- 
tained from elastic scattering of electrons from 
nuclei 7 . The charge distribution obtained from 
those experiments is often parametrized with the 
help of two-parameter Fermi model [3]: 



P ( r ) = Po 1 + exp 



(10) 



where c is the radius of the nucleus, a is the so- 
called diffiusness parameter of the charge density. 

Fig. [5] shows the charge density normalized to 
unity. The correct normalization is: po,A«(0) = 




Figure 2. The ratio of p the charge distibution to po 
the density in the center of nucleus. 
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Figure 3. The moduli of the charge form factor 
F em (q) of the 197 Au and 208 Pb nuclei for realistic 
charge distributions. For comparison we show the 
monopole form factor for the same nuclei. 



l -fm~ 



^f±fm~ 3 for Au and po iPb (0) = ^ 
for Pb. 

The form factor is the Fourier transform of the 
charge distribution [7]: 



f 47T 

F(q) = J — p (r) sin (qr) rdr 



(11) 
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Fig. [3] shows the moduli of the form factor as a 
function of momentum transfer. Here one can 
see many oscillations characteristic for relatively 
sharp edge of the nucleus. The results are de- 
picted for the gold (solid line) and lead (dashed 
line) nuclei for realistic charge distribution. For 
comparison we show the monopole form factor 
used in the literature. The two form factors co- 
incide only in a very limited range of q. 
The monopole form factor 



leads to simplification of many formulae for pro- 
duction of pairs of particles via photon-photon 
subrocess in nucleus-nucleus collisions. In our 
calculation A is adjusted to reproduce root mean 

perimental data |S]. 

3. Examples 

In the following we shall discuss several exam- 
ples considered by us recently We shall discuss 
individual processes in respective subsections. 

3.1. Exclusive production of pairs 

For dimuon production the elementary cross 
section can be calculated within QED. 

In Ref.[3] we have presented several distribu- 
tions in muon rapidity and transverse momentum 
for RHIC and LHC experiments, including ex- 
perimental acceptances. Here we wish to present 
only one example. 

The ALICE collaboration can measure only 
forward muons with psudorapidity 4 < rj < 5 
and has relatively low cut on muon transverse 
momentum p t > 2 GeV. In Fig|4] (left panel) we 
show invariant mass distribution of dimuons for 
monopole and realistic form factors including the 
cuts of the ALICE apparatus. The bigger invari- 
ant mass the bigger the difference between the 
two results. The same is true for distributions in 
muon transverse momenta (see the right panel). 

3.2. Exclusive production of 7r + 7r pairs 
Basic diagrams of the Brodsky and Lepage for- 
malism are shown in Fig. [5j The invariant ampli- 




ALICE conditions: 1/3,1/4 = (3,4), p3 t ,P4t > 2 GeV. 




Figure 5. Feynman diagrams describing the 77 
(qq)(qq) — > nn amplitude in the LO pQCD. 



tude for the initial helicities of two photons can 
be written as the following convolution: 



M(X 1 ,X 2 ) = 




(j>K (x, nl) T^ lA2 (x, y, fi 2 ) ^ (y, My) , (13) 

where \l x = min (x, 1 — x) y/s(l — z 2 ), fj, y = 
min (y, 1 — y) \/ s(l — z 2 ); z = cos# [12]. We take 
the helicity dependent hard scattering amplitudes 
from Ref. [13]. These scattering amplitudes are 
different for 7r + 7r _ and 7r°7r°. The distribution 
amplitudes are subjected to the ERBL pQCD 
evolution |15ll6j . The scale dependent quark dis- 
tribution amplitude of the pion can be expanded 
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in term of the Gegenbauer polynomials: 

oo 



(14) 



/„• above is the pion decay constant. 

Different distribution amplitudes have been 
used in the past. Wu and Huang [18] proposed 
recently a new distribution amplitude (based on 
a certain light-cone wave function): 



V3Am q l3 



2^3/2/, 
\ + Bx C*2 /2 (2a - 1) 



Erf 



Wxil-x) 



Erf 



8,3 2 :r (1 - x) 



(15) 



This pion distribution amplitude at the initial scale 
is controlled by the parameter B. They have found 
that the BABAR data for pion transition form factor 
at low and high transferred four-momentum squared 
regions can be described by setting B to be around 
0.6. This pion distribution amplitude is rather close 
to the well know Chernyak-Zhitnitsky [17] distribu- 
tion amplitude (4>ttcz — 3Cte(l — x)(2x — l) 2 ). In the 
following we shall use B — 0.6 and m q = 0.3 GeV. 
Then A = 16.62 GeV" 1 and j3 = 0.745 GeV. 

The total (angle integrated) cross section for the 
process can be expressed in terms of the amplitude of 
the process discussed above as: 



2tt 



P \ " 

4 ■ 64tt 2 W 2 q ^ 



|A4(Ai,A 2 )| 2 dz,(16) 



where the factor 4 is due to averaging over initial 
photon helicities. 

The hand-bag model was proposed as an alterna- 
tive for the leading term BL pQCD approach [19] . 
It is based on the philosophy that the present ener- 
gies are not sufficient for the dominance of the lead- 
ing pQCD terms. As in the case of BL pQCD the 
hand-bag approach applies at large Mandelstam vari- 
ables s ~ — t ~ — u i.e. at large momentum transfers. 
Diehl, Kroll and Vogt presented a sketchy derivation 
[19| obtaining that the angular dependence of the am- 
plitude is oc 1/ sin 2 9. 



In this approach the ratio of the cross section for 
the 7r°7r° process to the 7r + 7r~ process does not de- 
pend on 6 and is ^. The nonperturbative object 
R nw (s) in the hand-bag amplitude, describing tran- 
sition from a quark pair to a meson pair, cannot be 
calulated from first principles. In Ref. [19] the form 
factor was parametrized in terms of the valence and 
non- valence form factors as: 



(17) 



The a«, n u , a s and n s values found from the fit in 
Ref. 19 slightly depend on energy. For simplicity we 
have averaged these values and used in the present 
calculations: a u = 1.375 GeV 2 , n u — 0.4175, a s = 
0.5025 GeV 2 and n 3 = 1.195. 

In Ref. [6] we have discussed in detail elementary 
cross sections as a function of photon-photon energy 
and as a function of cos(6 l ). Here we have room 
for presenting only nuclear cross sections calculated 
within EPA discussed in the theoretical section. 

In Fig. [6] we show distribution in the two- pion in- 
variant mass which by the energy conservation is also 
the photon-photon subsystem energy. For this fig- 
ure we have taken experimental limitations usually 
used for the -kit production in e + e~ collisions. In the 
same figure we show our results for the 77 collisions 
extracted from the e + e~ collisions together with the 
corresponding nuclear cross sections for tv + tt~ (left 
panel) and (right panel) production. We show 

the results for the standard BL pQCD approach and 
for the hand-bag approach. 

By comparison of the elementary and nuclear cross 
sections we see a large enhancement of the order of 
10 4 which is somewhat less than Z\Z\ one could ex- 
pect from a naive counting. 

3.3. Exclusive production of p°p° pairs 

At low energies one observes a huge enhancement 
of the cross section for the elementary process 77 — > 
p°p° (see left panel of Fig[7]). In the right panel 
we show predictions of a simple Regge-VDM model 
with parameters adjusted to the world hadronic data. 
More details about our model can be found in our 
original paper [2]. 

In Fig[8] we show distribution in p°p° invariant 
mass (left panel) and the ratio of the cross section 
for realistic and monopole form factors. 

3.4. Exclusive production of cc 

In Fig l9ll0lllll2l we show photon-photon processes 
leading to the QQ in the final state. In the following 
we shall discuss them one by one. 
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Figure 6. The nuclear (upper lines) and elementary 
(lower lines) cross section as a function of photon- 
photon subsystem energy W- yl in the b-space EPA 
within the BL pQCD approach for the elementary 
cross section with Wu-Huang distribution amplitude. 
The angular ranges in the figure caption correspond 
to experimental cuts. 



Figure 7. Elementary cross section for 77 — > p°p° 
reaction. The fit to the experimental data is shown in 
the left panel and our predictions for the high energy 
in the right panel. 



Let us start with the Born direct contribution. The 
leading-order elementary cross section for 77 — > QQ 
as a function of W~ ll takes a simple form which dif- 
fers from that for 77 — > l + l~ by color factors and 
fractional charges of quarks. 

In the current calculation we take the following 
heavy quark masses: m c = 1.5 GeV, mi = 4.75 GeV. 
This formula can be directly used in the impact- 
parameter-space EPA. It is obvious that the final QQ 
state cannot be observed experimentally due to the 
quark confinement and rather heavy mesons have to 
be observed instead. Presence of additional few light 
mesons is rather natural. This forces one to include 
more complicated final states. 

In contrast to QED production of lepton pairs in 
photon-photon collisions, in the case of QQ produc- 
tion one needs to include also higher-order QCD pro- 
cesses which are known to be rather significant. Here 
we include leading-order corrections only for the di- 
rect contribution. In a s -order there are one-gluon 
bremsstrahlung diagrams (77 — >• QQg) and interfer- 
ences of the Born diagram with self-energy diagrams 
(in 77 — > QQ) and vertex-correction diagrams (in 
77 — ► QQ)- The relevant diagrams are shown in 
Fig llOl In the present analysis we follow the approach 
presented in Ref. [20]- The QCD corrections can be 
written as 



QCD 



(W 7 



N c e 4 Q 



~W^~ TV 



(1) 



(18) 
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Figure 8. Distribution in the p p invariant mass 



The function is calculated using a code provided 
by the authors of Ref. 20 . In the present analysis 
the scale of a s is fixed at p 2 = imq. 

We include also the subprocess 77 — > QQqq, where 
q (q) are u, d, s, quarks (antiquarks). The cross sec- 
tion for this mechanism can be easily calculated in the 
color dipole framework 21 22 . In the dipole-dipole 
approach [22] the total cross section for the 77 — > QQ 
production can be expressed as 



,0 (W, 



1,21. 



* /2/2 [P2,Z 2 ) 



Odd (pi, P2,XQf) d pidzid p2dz2 
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Figure 9. Representative diagrams for the Born am- 
plitudes. 
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Figure 10. Representative diagrams for the leading- 
order QCD corrections. 



E 



(P2,Z 2 ) 



<!-'.>' 

/- ./. /- ./.. (19) 



Udd{pi, P2,Xfo) d p\dz\d p2dz 2 , 



where (p, z) are the quark - antiquark wave 

functions of the photon in the mixed representation 
and add is the dipole-dipole cross section. Eq.(|19[) 
is correct at sufficiently high energy W 77 3> 2mQ. 
At lower energies, the proximity of the kinematical 
threshold is a concern. In Ref. [5T] a phenomenologi- 
cal saturation-model inspired parametrization for the 
azimuthal angle averaged dipole-dipole cross section 
has been proposed: 



a dd = cr o 



1 



exp 



4i?g (x a 



Here, the saturation radius is defined as 

-A/2 



Ro {Xab) 



1 (Xab\ ' 

Qo V X ) 



(20) 



(21) 




Figure 11. Representative diagrams for QQqq pro- 
duction. The oval in the figure means a complicated 
interaction which is described here in the saturation 
model as explained in the main text. 



and the parameter x a b which controls the energy de- 
pendence was given by 



Xab 



4m„ + Am 2 b 



(22) 



The effective radius is parametrized as r cff — 
(P1P2) 2 /(pi + P2) [21] . Some other parametrizations 
of the dipole-dipole cross section were discussed in 
the literature. The cross section for the 77 — > QQqq 
process here is much bigger than the one correspond- 
ing to the tree-level Feynman diagram as it effectively 
resums higher-order QCD contributions. 

As discussed in Ref. 22 the QQqq component have 
very small overlap with the single-resolved compo- 
nent because of quite different final state, so adding 
them together does not lead to double counting. The 
cross section for the single-resolved contribution can 
be written as: 



dxi [gi ( 



Xl,fl 



XiS 



)] + 



dx 



2 52 



(x2,p 2 ) <j lg (I = X2S)] , 



(23) 



where g\ and g2 are gluon distributions in photon 1 
or photon 2 and <r 97 and o lg are elementary cross 
sections. In our evaluation we take the gluon distri- 
bution from Ref. 1231. 
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Figure 12. Representative diagrams for the single- 
resolved mechanism. The shaded oval means either 
t- or u- diagrams shown in Fig. |9] 



Elementary cross sections have been presented and 
discussed in Ref. [4]. Here we show only nuclear cross 
sections. In Fig. [13] we compare the contributions of 
the different mechanisms as a function of the photon- 
photon subsystem energy. For the Born case it is 
identical as a distribution in quark-antiquark invari- 
ant mass. In the other cases the photon-photon sub- 
system energy is clearly different than the QQ in- 
variant mass. These distributions reflect the energy 
dependence of the elementary cross sections. Please 
note a sizable contribution of the leading-order cor- 
rections close to the threshold and at large energies 
for the cc case. Since in this case VFy 7 > Mqq, 
it becomes clear that the QQqq contributions must 
have much steeper dependence on the QQ invariant 
mass than the direct one which means that large QQ 
invariant masses are produced mostly in the direct 
process. In contrast, small invariant masses (close 
to the threshold) are populated dominantly by the 
four-quark contribution. Therefore, measuring the 
invariant mass distribution one can disentangle some 
of the different mechanisms. As far as this is clear for 
the cc it is less transparent and more complicated for 
the bb production. In the last case the experimental 
decomposition may be in practice not possible. 

Finally in Table 1 we show partial contribution of 
different subprocesses discussed above. 

3.5. Final comment 

We have presented four examples of processes that 
could be studied at RHIC or LHC. In all cases size- 
able cross sections have been obtained. Correspond- 
ing measurements are not easy as one has to assure 
exclusivity of the process, i.e., it must be checked that 
there are no other particles than that measured in 
central detectors. In all cases fissibilty studies, in- 




20 40 60 80 100 20 40 60 80 100 



W„ (GeV) W„ (GeV) 



Figure 13. The nuclear cross section as a function of 
photon-photon subsystem energy W~ ll in EPA. The 
solid line denotes the results corresponding to the 
Born amplitude (cc -left panel and bb -right panel). 
The leading-order QCD corrections are shown by the 
dash-dotted line. For comparison we show the differ- 
ential distributions in the case when an additional 
pair of light quarks is produced in the final state 
(dashed lines) and for the single-resolved components 
(dotted line). 



Table 1 

Partial contributions of different mechanisms at 



\J snn — 5.5 TeV. 





Ctot 


Born 


QCD-corr. 


4-q 


Sin.-res. 


cc 


2.47 mb 


42.5 % 


14.6 % 


27.1 % 


15.8 % 


bb 


10.83 [ib 


18.9 % 


7.7% 


64.5 % 


8.9 % 



eluding Monte Carlo simulations, are required. 
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